- 1 - 



Numerical Analysis of Magnetic Field Amplification by 

Turbulence 

by 

Hongsong Chou 

Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138, U.S.A. 

chou5@f as. harvard, edu 

Received ; accepted 



- 2 - 
ABSTRACT 

We apply a Fourier spectral numerical method to 3D incompressible MHD 
turbulence with a magnetic Prandtl number Pr > 1. We examine the processes 
by which an initially weak, large-scale seed magnetic field and an initially 
weak, small-scale, impulse-like seed magnetic field are amplified. We find that 
in both cases the magnetic energy spectrum grows at all scales. The growth 
rates at different amplification stages are analyzed. For a large-scale seed 
magnetic field, the magnetic energy density grows as ~ t 2 for the first few 
turbulence eddy turnover times, followed by a dynamic growth stage, where 
nonlinear interactions between different scales of the turbulence contribute to 
an exponential growth rate that is largely determined by the turbulence eddy 
turnover time. For a seed magnetic field that is initially set up at a small scale 
in the turbulence, during the kinematic development stage, the growth rate 
of magnetic energy is oc l/r max , where r max is the eddy turnover time of the 
smallest eddies of the turbulence. The kinematic growth stage is followed by a 
dynamic growth stage, where nonlinearity plays important role. During such 
dynamic growth stage, the growth rate of total magnetic energy is determined 
by both the magnetic energy amplification within the turbulence inertial range 
and that within the turbulence dissipation range. 

Subject headings: MHD turbulence, Dynamo 
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1. Introduction 

Astrophysical magnetic field has often been observed on much larger scales than the 
scales of astrophysical turbulence. The connection between the generation of large-scale 
magnetic field and the turbulence of much smaller scales has been contemplated by 
astrophysicists for many decades. Mean-field electrodynamics (MFE, see Moffatt 1978 
or Krause & Radler 1980), among other theories, employs a two-scale approach to the 
problem. It suggests that the helical turbulent motions of astrophysical plasma may align 
small-scale magnetic field so that an observable, large-scale magnetic field can be formed. 
Because MFE is largely a linear theory, its applicability to nonlinear phenomena such as 
the solar or interstellar turbulence has been questioned ever since the introduction of this 
theory. Some early criticism of MFE was discussed by Piddington (1975). He argued that 
kinematic solar dynamo theories do not account for the removal of the large amounts of 
flux generated each solar cycle. Recent objections to dynamo action have their root in the 
problem of small-scale magnetic fields. The amplification of seed magnetic field in galaxies 
has been considered by Kulsrud & Anderson (1992). The magnetic Prandtl number, defined 
as Pr — vj\ where v is the molecular viscosity and A the magnetic resistivity, is much 
greater than 1 in the interstellar medium of galaxies. Kulsrud & Anderson predict that 
the growth rate of an initially weak, small-scale magnetic field in Pr ^> 1 systems will 
be oc l/r max , where r max being the eddy turnover time of the smallest turbulent eddies. 
Because r max is very small in astrophysical turbulence, the growth rate of the small-scale 
magnetic field will be large. They then argue that because of such fast amplification of 
magnetic energy at small scales, the turbulence may be reduced dramatically in a short 
period of time, so that it is unable to amplify the magnetic field at scales much larger than 
the turbulent dissipation scale. The recent work by Kinney et al. (2000, see also Cowley 
2000) also criticizes the galactic dynamo theory by arguing that because the galactic plasma 
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has a very large magnetic Prandtl number, any small-scale seed magnetic field will grow 
quickly to lock the velocity field in a shear motion pattern, where the dissipation term in 
the momentum equation, z/V 2 U, is balanced by the Lorentz force term, B • VB. Here U 
and B are the velocity field and the magnetic field. 

This work is motivated by the studies of Kulsrud & Anderson and Kinney et al. 
We noticed that the work by Kulsrud & Anderson can be valid only for the kinematic 
development of the magnetic field, as back reaction of the magnetic field on the velocity 
field was not considered in their governing equations. Therefore, their prediction that the 
magnetic field at large scales may not grow after the velocity field is quenched at small 
scales needs examination. We also noticed that the simulation by Kinney et al. is for 2D 
MHD, which can be very different from the case of 3D MHD. Moreover, they focused their 
research on those scales that are smaller than the velocity dissipation scales, and did not 
include the inertial range of the turbulence. To further understand the physics of magnetic 
field amplification by turbulence and its application to 3D astrophysical systems, we carried 
out a numerical study of incompressible 3D MHD systems. With our numerical model, we 
study the amplification of initially weak, large-scale and small-scale seed magnetic field by 
turbulence. Unlike the theoretical analysis of Kulsrud & Anderson, our numerical study is 
not restricted to the kinematic development stage of the MHD system. Rather, we can also 
study the effects of back reaction. Unlike the numerical work by Kinney et al. (2000), we 
include an inertial range in our numerical model, and study whether or not the magnetic 
field within the inertial range will grow, especially when the magnetic field beyond the 
turbulence dissipation scale {I v ,d) grows and significantly modifies the velocity field near 

lv,D- 

The structure of this paper is as follows: we introduce our numerical model in section 
2; sections 3 and 4 are devoted to the detailed numerical analysis of magnetic energy 
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spectrum development for a large-scale seed field and a small-scale seed field, respectively; 
the physical implications and applications of our numerical results are discussed in section 
5; conclusions and suggestions for future work are given in section 6. Our numerical study 
is different from those of Brandenburg (2000) and Cho & Vishniac (2000) in the following 
ways: first, not like the work of Cho & Vishniac, in most of our simulations, we have a 
magnetic Prandtl number Pr = 3, so that the velocity dissipation scale is larger than the 
magnetic dissipation scale; second, there is no scale separation as the one discussed in 
Brandenburg (2000); third, although Brandenburg (2000) studied cases of various magnetic 
Prandtl numbers, the initial conditions for the simulation runs in that work are different 
from those in our work; therefore, the physical processes considered in that work and ours 
are different. 

In this work, we distinguish four scales. We denote the size of turbulence energy 
containing eddies as L. L is also called the outer scale or the integral scale of the turbulence. 
We denote the dissipation scales of the velocity field and the magnetic field as I v> d and I^d, 
respectively. For turbulence of large kinetic Reynolds number, we have L ^> l Vt o- In Fourier 
space, we introduce two wave numbers that correspond to I v> d and I^d' K,d ~ ^-/K,d and 
kb,D ~ l/h,D- In some of our simulation runs, we have kb,D > k v ,D, i-e., the magnetic 
dissipation scale is smaller than the velocity dissipation scale. This is consistent with a 
magnetic Prandtl that is greater than 1. In other simulation runs, we have Pr = 1. Another 
scale is the so-called "ensemble average scale" , which is denoted by T ^> L and over which 
we calculate averaged quantities, < • > (or ~). Finally, we denote the scale of the whole 
physical system as S, which is the typical scale for the variations of averaged quantities 
< • > (or ~). What we mean by "large scale" in the following sections is the scale ~ S. So 
we have the relation S » T » i > l V:D > l b>D . In our model, there are two large-scale 
quantities, B and V. We assume both of these quantities are constant, i.e., S — > oo. We 
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set our reference frame to that moving at V and henceforth omit terms of V. To aid our 
discussion, we decompose the total magnetic field into two parts: B = B + b, where b is 
the fluctuating component of B. To make our discussion easy to follow, throughout the 
paper we will use a few terms and notations, which will be introduced as we proceed. They 
are also listed in Table I. Readers may refer to this table for further clarification. 



2. The Numerical Model 

We focus on the incompressible MHD equations, which are solved numerically using the 
standard Fourier spectral method. The nonlinear terms are evaluated by a pseudo-spectral 
procedure. Let B and U be the magnetic field and the velocity field, respectively. Under 
an external forcing term J 7 , the undimensionalized incompressible MHD equations can be 
written as (with Einstein summation convention) 

(d t - iCV 2 ) Ut = d 3 {-p5ij - U t U 3 + BiB 3 ) + T h (1) 

(d t - iCV 2 ) B t = d 3 (UiBj - BiUj) , (2) 

diUi = d 3 B 3 = 0. (3) 

where R e and R m are the kinematic and magnetic Reynolds numbers and defined as 

R e = V -^,R m = V -^. (4) 
v A 

Here v rms is the root-mean-square of the velocity field, L the integral scale of the turbulence, 
v the molecular viscosity, A the magnetic resistivity. Note that we have already written B 



in units of y/Airp after we divide both sides of the momentum equation and the induction 
equation by density p. If we use a hat, A, to denote discrete Fourier transform, and <S> to 
denote convolution, the above equations in Fourier space are 

fa + vk 2 ) Uj = Pji [ik m (-U t (g)U m + B l (g) B m ) + JFj , (5) 



(d t + AA; 2 ) Bj = Pji [ik m (U t <g> B m -U m ®Bi)], (6) 

k m U m = k j B j = 0. (7) 

Here P is the projection operator defined as Pji = 5ji — In our simulation, we treat the 
system as a cube [0, 2n) x [0, 2n) x [0, 2ir). The Cartesian coordinate of a grid point can be 
written as xi = jf-l, y m = j^m, z n = jf-n, for l,m,n — 0, 1, 2, JV — 1. A point in Fourier 
space has coordinates k s = s,k p = p, k q = q, for s,p,q— — y, — y — 1. 

Equations (5), (6) and (7) are numerically solved with the standard Fourier pseudo- 
spectral method. Equations (5) and (6) are treated as ordinary differential equations for U 
and B. With the projection operator P, the divergence free condition (7) will be satisfied 
for t > as long as U and B are divergence free at t — 0. All our simulations start from 
divergence free initial conditions. We employ a second-order Runge-Kutta (RK2) method 
to advance equations (5) and (6) in time. We can exploit the advantage of using RK2 in the 
following two aspects. First, an integral factor can be easily introduced with the transform 

U m (t) = U m (K t)e~ uk2 \ B m (t) = S m (k, t)e- Xk2t . (8) 

Second, aliasing errors can be reduced by introducing positive and negative random phase 
shifts at the first and second stages of RK2, respectively (Machiels & Deville, 1998). The 
forcing term used in our simulation is a combination of the one used by Chen et al. (1993a) 
and the one used by Brandenburg (2000), with slight modification. Chen et al. employed 
a forcing term that maintains the energy density values at k — 1 and k = 2 so that the 
energy spectrum of velocity field always follows Kolmogoroff k~ 5 ^ 3 law at k — 1 and k = 2. 
In Brandenburg's simulation, he injected helical waves of random phases at k — 5. Because 
our simulations start from an established hydrodynamic turbulence, we first use the forcing 
by Chen et al. to obtain a fully developed pure hydrodynamic turbulence. Such fully 
developed turbulent velocity field is taken as the initial velocity field for the following MHD 
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turbulence simulation. For all the MHD turbulence simulation runs, the forcing function 
has the form 



That is, the force works only within the shell S\ : 0.5 < |k| < 1.5. Here T c is a forcing 
term that is similar to the one adopted by Chen et al. It is calculated by multiplying 
the velocity components within shell S\ by a factor, 7 > 1, so that before a new step of 
integration starts, the kinetic energy density within this shell is reset to E\ = 0.24. Phases 
of the velocity components within the shell are not changed. This forcing is equivalent as 
lengthening the velocity vector within shell Si by a factor 7 — 1. Denote the increment of a 
velocity vector under force T c as 5v = R + il, where R, I are the real and imaginary parts 
of SU. In a few runs of our simulation, we need to inject kinetic helicity into the turbulence. 
To do this with JF C , we tune the angle between R and I so that they remain perpendicular 
to each other. Because kinetic helicity at k can be calculated as H(k) = 2k • R x I, in doing 
so, we inject kinetic helicity into the turbulence. 

The forcing term T c maintains the energy level at the forcing scale so that the 
fluctuation in the energy development history can be small; therefore, the growth stages of 
both the kinetic energy and the magnetic energy can be studied carefully and accurately. 
However, this force does not introduce random phases into the velocity field. To be more 
realistic about the forcing in our simulation, we also use the forcing term derived 
from the forcing function used by Brandenburg (2000), as a secondary forcing function to 
introduce random phases into the velocity field. JF 6 has the form 



Here JF < 1 is a factor adjusted at each time step so that the kinetic energy density 
within shell Si fluctuates within 5% of E\. e is an arbitrary unit vector in Fourier space. 



^(0.5 < |k| < 1.5) = T c + T h . 



(9) 




(10) 
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4>(t) is a random phase. Note that jF(k)* = JF(— k) so it is real, and it is helical in that 
J= ■ V x T = -kT 2 < 0, i.e., it has maximum helicity. Because Th is tuned in such way that 
it only contributes to 5% of the kinetic energy at the forcing scale, can be considered as 
a perturbation to JF C . Therefore, the advantage of using (9) as the forcing term is three fold: 
to avoid strong fluctuations of kinetic and magnetic energy density with time, to introduce 
random phases to the velocity field, and to maintain the kinetic helicity at certain level. 

In order to study the non-unit Prandtl number case, we adopt the following hyper- 
viscosity and hyper-diffusivity: — z/ 7 /c 14 U and — A 7 /c 14 B. The dissipation scale of the velocity 



of the hyper-viscosity in our simulations. e v is the mean dissipation rate of the velocity 
field and is calculated using the formula e v = 2v h J2kk 2h E v (k), with h = 7. E v (k) is the 
kinetic energy spectrum, k^p is defined in a very similar way as the one k v> r> is defined 
above. With z/ 7 = 5.0 x 10~ 16 and A 7 = 3 x 10~ 20 we have the dissipation scales of velocity 
field and magnetic field as k v> n ~ 13, k^o ~ 28, respectively, so that the Prandtl number is 
Pr = (^y 6 » 3. In many other studies (see Brandenburg 2000), the magnetic Prandtl 
number is defined as the ratio of molecular viscosity to magnetic resistivity. Because we 
apply hyper-viscosity in most of our simulations, the magnetic Prandtl number used in this 
work can be considered as an "effective" magnetic Prandtl number, i.e., a ratio inferred 
from the measured dissipative cutoff wavenumbers. 

A ID version of our code is used to solve the nonlinear Burger's equation (see section 
6.1 of Canuto et al. 1988) and the numerical results match exactly the analytic results. If 
we impose a strong uniform magnetic field and a small disturbance of velocity field, a pair 
of Alfven waves are numerically generated, both propagating along the uniform magnetic 
field but in opposite directions. This is exactly predicted by linearized incompressible MHD 
equations. Finite amplitude MHD waves are also produced in our numerical simulations (see 




(Machiels & Deville, 1998). Here h = 7 is the order 
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section 10.1 of Moffatt, 1978). Our code is also used to study hydrodynamic turbulence with 
normal dissipation. With the forcing given by (9), the flows reach a statistically stationary 
state in five to ten large eddy turnover times. In Figure 1, we plot the kinetic energy spectra 
of two different pure hydrodynamic turbulence simulation runs with normal dissipation. 
The kinetic energy spectrum is calculated as E v (k) = \ J2k-ol v(k') 2 . The magnetic 
energy spectrum shown in next few sections is calculated as E b (k) = \ YX-q\ H^') 2 - I n 
Figure 1, R\ — v rms A/u is the Taylor micro-scale Reynolds number. A is the Taylor 
micro-scale defined by A = ^Jl5uv 2 ms /e, where e is the rate of dissipation of hydrodynamic 
kinetic energy per unit mass (Chen et al. 1993b). For R\ = 43, the exponential falloff 
starts around k ~ 6.5, while for R A = 70, the falloff starts around k ~ 11. Both cases 
exhibit a Kolmogoroff k~ 5 ^ 3 inertial range. With a hyper-viscosity, the turbulence in 
stationary state will have a energy spectrum that deviates from a Kolmogoroff k~ 5 ^ 3 inertial 
range. Rather, the spectrum will be flatter than k~ 5 ^ 3 , as shown in numerous simulations 
(Michiels & Deville, 1998). We plot the kinetic energy spectrum with the hyper- viscosity 
u 7 = 5.0 x 10~ 16 in Figure 1. The spectrum in the inertial range follows a power law of 
~ AT 12 . In Table 2, we have listed the parameters of all the simulation runs that we have 
done for this work. The resolution of our simulations is (64) 3 . 



3. Amplification of a large-scale seed magnetic field by the turbulence 

The amplification of an initially weak, large-scale magnetic field is studied by imposing 
a constant magnetic field along y-direction, B = By with B = 0.0316, into a homogeneous, 
isotropic and stationary hydrodynamic turbulence under the forcing of (9). The magnetic 
energy density associated with this initial field is Eb = 5 x 10~ 4 . Because the turbulence 
will stretch B, magnetic field at k > 1 will be generated. If the magnetic field is weak, 
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its back reaction on the velocity field is small; therefore, the turbulence maintains its 
stationarity until the magnetic field grows strong enough to alter the flows. The growth 
rate of magnetic field due to the presence of a constant B with negligible back reaction can 
be estimated as follows. The induction equation can be written as 

d t b = (B • V)U + AV 2 b. (11) 

Therefore, shortly after the start of the simulation, i.e., during the first few eddy turnover 
times, the magnetic energy spectrum within the inertial range can be calculated as 

|b(k)| 2 ^t 2 5 2 A; 2 |U(k)| 2 . (12) 

If the kinetic energy spectrum follows |U(/c)| 2 oc k~ p , we have |b(k)| 2 oc k~ p+2 . Figure 2 
shows the energy spectra of the velocity field and the magnetic field at 0.1r e ddy after the 
start of simulation. The kinetic energy spectrum changes little, while the magnetic energy 
spectrum in the inertial range follows a power law of A; 77 , as predicted. To show how fast 
the magnetic energy grows in this stage, we calculate the following quantity 

m = ~ - (is) 

1 y ! At E b t K ' 

for all fc's as well as the total magnetic energy. f3(t) can be considered as an instantaneous 
growth rate. It is an explicit function of and varies with time. From equation (13), we find 
that (3{t) should not be a function of k, and this is clearly shown in the left panel of Figure 
3, where we plot f3(t) at different fc's. At different fc's, (3(t) is roughly the same, and decrease 
with time. Near dissipation scales, (3{t) decreases with time slightly faster than those within 
the inertial range, and this is due to the strong dissipation at large fc's, i.e., the dissipation 
term in (11) becomes important. In the right panel of Figure 3, we plot the averaged value 
(over k) of f3(t), and it behaves like | shortly after the start of simulation, but decreases 
faster than f as the strength of b grows to be comparable with the strength of B. This 
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confirms equation (13). The deviation of the averaged f3(t) from | after t = 4 becomes more 
and more prominent as the nonlinear interaction between the velocity field and the growing 
magnetic field becomes stronger and stronger. This is shown in Figures 4 and 5. From 
Figure 4, we find that the growth of the magnetic energy density follows four stages. Stage 
1 starts from the beginning of the simulation till t = 3, during which the magnetic energy 
density grows as t 2 . This is a linear stage in that the line stretching of B by the velocity 
field contributes to most of the growth of the magnetic energy, while the velocity field 
changes little. During stage 2 that lasts from t = 4 till t = 16, it grows exponentially with a 
growth rate (5 — 0.1. The velocity field is suppressed dramatically by the end of stage 1 and 
at the beginning of stage 2. Figure 4 shows that the kinetic energy density drops by 25% 
during this period. After the growth of magnetic energy enters the exponential stage, the 
loss of kinetic energy slows down. Within stage 2, the velocity field loses 5% of its initial 
energy at t — 0. The slowdown of the energy loss of velocity field is due to the existence 
of a forcing, which injects energy at 0.5 < |k| < 1.5 so that the £"^(0.5 < |k| < 1.5) is 
maintained at 0.240 ± 0.012. Because the nonlinear interaction between the velocity field 
and the magnetic field is strongest near the velocity dissipation scale k v ^o, the velocity 
at small scale is suppressed most, which can be seen from Figure 5. The kinetic energy 
density at k v ^ D = 13 drops an order of magnitude from t = 0.31 to t = 9.42. Therefore, 
the continuing growth of the total magnetic energy during stage 2 should be attributed to 
the energy input from the forcing scale and the velocity line stretching within the inertial 
range. During stage 2, the dominant terms in the induction equation within the inertial 
range is the line stretching term b • VU, and this is because the fluctuating component b of 
the total magnetic field has grown greater than the imposed B. The nonlinear interaction 
between U and b, and between b at different scales provides a self-excitation of b that 
is independent of B and is capable of exponentially amplifying b. Such an exponential 
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growth rate should be determined by the statistical properties of the velocity field. In our 
simulation, the largest integral length that the largest eddy circles around can be calculated 
as L = 2nJ2kk~ 1 E v (k)/J2kE v (k). The rms velocity of the turbulence can be estimated as 
v rms = \J^J2k E v (k). Within stage 2, we found that the temporal average of L (denoted 
by < L > t ) is < L > t ~ 4.62, and the temporal average of v rms (denoted by < v rms > t ) is 
< v rms > t ?s 0.52. Therefore, l/r e ddy =< v rms > t / < L > t = 0.11 ~ (3, which is consistent 
with the theoretical prediction of Parker (1979). 

After the exponential growth stage, magnetic energy density growth enters a near 
saturation stage from t — 16 to t — 30, where the growth is further slowed down. During 
this stage, the magnetic energy at small scales almost stops growing, while at large scales 
the magnetic field continues to grow. For the saturation stage, magnetic energy density 
fluctuates around 0.33, while the kinetic energy density fluctuates around 0.4. 

The growth of magnetic structures can be seen from Figure 6, where we plotted 
isosurfaces of the magnetic field strength, |B| = 2.5(|B|). Here (|B|) is the spatial mean 
of the magnitude of the magnetic field. The isosurface at a time shortly after the start of 
the simulation is shown in panel (a) of Figure 6. At this time point, the magnetic field 
consists of mostly small-scale structures. The structures shown in panel (a) occupy 3.5% of 
the whole system and are distributed quite evenly in space. However, at a later stage of the 
development, the magnetic field develops structures that span scales that are close to the 
forcing scale, k as 1, in the system. This is shown in panel (b) at t — 31.4 of the saturation 
stage. Tube-like structures that span the whole simulation box can be clearly seen. 

In our simulation of Run A, the kinetic helicity, calculated as < v ■ V x v > as 
a function of time, has a value of H v = —0.13 ± 0.04. Most of the kinetic helicity 
comes from the forcing scale, as the force is helical. The normalized kinetic helicity is 
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H v / (tu rms ■ v rms ) « -0.13/0.74 • 0.51 = 0.27 ± 0.08. Such helical flow will generate a 
dynamo a— effect and the a coefficient can be estimated as a ~ —i^- < v ■ V x v >~ 0.05, 
where r cor is the correlation time of velocity field and is measured to be r cor ~ 1.0. A growth 
rate of the magnetic energy density due to such a dynamo a— effect can be estimated as 
P a ~ 2a/ L ~ 0.02, which is smaller than the growth rate in the second stage, P = 0.1. 
Therefore, it is not clearly evident from this study that the dynamo a— effect in the 
moderately helical flow of our simulation is capable of driving the exponential growth of 
magnetic field in stage 2. Cho and Vishniac (2000) have done recent numerical simulations 
and also claim that a dynamo a— effect may play much less important roles in amplifying 
magnetic field in turbulent flows than the turbulent line stretching effect. One should also 
notice that the dynamo theory of mean-field electrodynamics requires scale separation 
between the outer scale of the turbulence and the scale of the physical system, while in our 
simulation there is no scale separation; therefore, one should use caution when applying 
our simulation results to the two-scale discussion of the dynamo a— effect. For two-scale 
approach to the dynamo problem, the reader is referred to a recent work by Brandenburg 
(2000), which presents a detailed account of large-scale magnetic field amplification through 
dynamo a— and P~ effects. 



4. Amplification of a small-scale seed magnetic field by the turbulence 

In this section, we present the numerical analysis of a seed magnetic field that is 
initially concentrated at a small scale, k = 20. We start from a fully developed turbulence 
and a magnetic impulse 

e for k = 20, 

E b {k) = { (14) 
otherwise. 
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Such an initial condition is consistent with the theoretical analysis by Kulsrud & Anderson 
(1992) on the amplification of weak small-scale magnetic energy, which they called magnetic 
noise. The growth of magnetic spectrum with eo = 0.001 (Run B) is shown in Figure 7. 
The initial impulse-like small-scale seed magnetic field is located at k = 20. After the 
simulation starts, the narrow impulse quickly becomes broader and broader, extending 
to both larger and smaller scales than k = 20. It extends to large wave numbers (small 
scales) and soon hits the magnetic dissipation scale k^D = 28 and the energy is removed 
by dissipation near k b) D- The broadened, impulse-like, seed field also extends to small 
wave numbers, transporting magnetic energy to larger and larger scales. From t = 2.2 
to t = 19.0, the portion of the growing magnetic spectrum that extends from k — 1 to 
k v> D follows a A; 3 / 2 power law, indicated by a group of dashed lines in Figure 7. Such k 3 ! 2 
profiles are reminiscent of the kinematic analysis of the magnetic spectrum by Kulsrud & 
Anderson (1992). Note that by Kulsrud & Anderson, such k 3 ! 2 profiles should extend into 
the range beyond k v> D but before k b ,D is reached. However, because of the low resolution of 
our simulation, we do not have a long range between k v D and k^D] therefore, the magnetic 
energy spectrum beyond h v ,D can be strongly affected by the dissipation process near the 
magnetic dissipation scale and does not exhibit the k 3 ! 2 scaling. Future simulations with 
higher resolution should be able to resolve such scaling within that range. 

To study the broadening of the initial magnetic impulse, we define the following two 
terms 



C b (k) 



E 



B*(p) 



Q(p), 



(15) 



fc-0.5<|p|<fc+0.5 



L b (k) 



E 



B*(P) 



R(p). 



(16) 



fc-0.5<|p|<fc+0.5 



Here 



Q = -U VB 



(17) 



16 



is the convection term, and — B ■ Q is the energy per unit time that is transferred away 
from magnetic field at one location to other locations. Also, 



is the term representing the line stretching effect, and B • R is the work done to the 
magnetic field by the velocity field. In Figure 8, we plot kC b {k) and kL b {k) at different time 
points. At t — 0.2, panel (a) of Figure 8 shows that the magnetic energy of the impulse at 
k = 20 is being transported to scales larger and smaller than k = 20. In panel (b) of Figure 
8, the energy is being removed from k = 20 to larger and smaller scales. It shows that the 
line stretching effect is generating magnetic structures at other scales. Panels (a) and (b) 
together explain the broadening of the impulse shown in Figure 7. Notice that at t = 0.2, 
the convection effect, i.e., C b (k), is comparable to the line stretching effect, i.e., L b {k). 
However, from panels of the right column of Figure 8, we find that the line stretching effect 
soon becomes dominant over the convection effect, which is shown at different times in 
the left column of Figure 8. kL b (k) is non-negative at all /c's, therefore always converting 
kinetic energy to magnetic energy. The convection term at different A;'s has different signs 
at different times, meaning that the magnetic energy can be transported into or out of 
certain scale of the MHD turbulence. Panels (b), (d), (f) and (h) also show that the line 
stretching effect is most prominent at small scales, that is, the scales between the velocity 
dissipation scale, k V; c = 13, and the magnetic dissipation scale, k bt r> = 28. 

According to Kulsrud & Anderson(1992), before the back reaction of magnetic field is 
strong enough to alter the velocity field, the exponential growth of magnetic energy from 
a magnetic impulse should be attributed to the line stretching at the velocity dissipation 
scale, i.e., 



R = B VU 



(18) 




(19) 
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where 7 = l/r max is the inverse of the eddy turnover time of velocity field at dissipation 
scale k V; o- Recent work by Chandran(1997) and Schekochihin & Kulsrud(2000) modified 7 
from 1/T max to e/r max and e « 0.6 is due to the non-zero correlation time of the velocity 
field. We have found such kinematic exponential growth stage in our simulations, such as 
the e^ 1 * portion of the magnetic energy growth in Figure 9. To test if the growth rate, f3±, 
in such exponential growth stage can be estimated according to the theory of Kulsrud & 
Anderson, we have measured this quantity from our simulation results. In Figure 10, we 
plot the kinematic exponential growth stage of the magnetic field, which extends from t — 2 
to t — 9 for Run C with an e (k = 20) = 10" 5 . The best fitting to E b (t) gives a growth 
rate of 7^ = 0.36 ± 0.05. We also calculated the growth rate due to line stretching by the 
largest eddies, ^eddy = 0.18 ± 0.04. The smallest eddy turnover time calculated from our 
numerical results within this time range is 7 n = 0.69 ± 0.12. Therefore, in our simulation, 
the modification factor due to non-zero correlation time is e = ^b/ln ~ 0.52, and such result 
is consistent with the work by Chandran (1997) and that of Schekochihin & Kulsrud (2000). 

For a Kolmogoroff turbulence, the eddy turnover time at scale k v> D can be calculated as 



With T e ddy an d k e ddy ~ 1/L from our simulation, we have ^Koi,D = ^/ T Koi,D = 0.57 < 
7„ = 0.69. This is due to the fact that the hyper-viscosity flattens the kinetic energy 
spectrum, hence produces larger values of velocity near the dissipation scale than that in 
Kolmogoroff turbulence. However, the nature of the growing of the magnetic field due to 
line stretching near the dissipation scale does not change with the slightly flattened kinetic 
energy spectrum. This can be seen from the work by Kulsrud & Anderson (1992). The 
growth rate of magnetic energy depends on the integral 




(20) 




(21) 
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where J(k, 0) kinetic energy density at zero frequency. For the velocity field in our 
simulation, we have J(k, 0) ~ k~ 1,2 , thus most of the contribution to 7 still comes from the 
smallest eddies. 

5. Discussion 

Our simulations show that the interaction between fully developed, constantly forced 
turbulence and an initially weak seed magnetic field will always lead to the growth of 
magnetic field at different scales of the turbulence. Given an initially weak, large-scale 
external magnetic field, the emergence of magnetic energy at small scales is due to the line 
stretching of B by the velocity field. This is a kinematic process, as the velocity field is not 
affected much by the growing magnetic field. Within the inertial range, where dissipation 
is negligible, the magnetic field is amplified mainly by the line stretching term, i.e., the 
second term of equation (11). This kinematic stage finishes when the nonlinear interaction 
between the growing magnetic field becomes comparable to the line stretching term, and a 
dynamic growth stage follows. In other words, when the terms U • Vb and b • VU are of 
the same order as B • VU, (11) is not valid anymore, and the self-excitation of the growing 
magnetic field within the inertial range will dominate other amplification processes. This is 
shown clearly as stage 2 in Figure 4, where the total magnetic energy grows exponentially. 
And from Figure 5, we find that such an exponential self-excitation stage happens after the 
magnetic energy densities at different scales grow to be comparable to the magnetic energy 
density of B. 

Several authors have been arguing that there might be a relation between the strength 
of a large-scale magnetic field, B, and the strength of the fluctuating component of the 
magnetic field, b. Krause and Radler (1980, chapter 7) derived a relation for 3D MHD 
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between B and (b 2 ) in the form 

(b 2 ) = B 2 ^, (22) 

where tjt and i] are the turbulent diffusion coefficient and the magnetic diffusivity, 
respectively. The governing equation they used in their derivation of (22) was equation (11), 
and the nonlinear terms Q, R in (17) and (18) were ignored. Our simulations show that 
this is not a valid procedure, as terms like Q and R are of importance to the exponential 
growth of the magnetic energy. In fact, from our simulations we find that no matter how 
large or small the initial large-scale magnetic field is, after 5 — 10 eddy turnover times of 
the turbulence, a statistically stationary magnetic energy spectrum will be formed. 

The amplification of the small-scale seed field, simulated in Runs B and C, further 
shows that relation (22) is not necessarily true, as B = in these two runs. The growth 
of magnetic field at all scales (see Figure 7) can be due only to the nonlinear interactions 
between the velocity field and the magnetic field at different scales. 

In the two simulation runs (B & C) with a moderate magnetic Prandtl number 
(Pr = 3), the magnetic energy within the inertial range, Ej, which is initially set to zero, 
grows to a steady state in which Ej > E v , where E v is the magnetic energy stored between 
the two dissipation scales, k Vj D and kb,D (see Table 1 for more details). Some authors 
(Cowley 2000) argue that for large magnetic Prandtl number, the magnetic field at small 
scales, i.e., k Vt o < k < k^o, grows so fast and strong that it swamps the velocity field at all 
scales and the magnetic field within the inertial range cannot grow at all. Our simulation, 
which includes both the inertial range and the dissipation range, does not support such 
picture. Instead, our simulation shows that Ej does grow. The initial magnetic impulse at 
small scales spreads to all scales between k Vt o and kj,^. The growth rate of magnetic energy 
at these small scales is approximately the inverse of the eddy turnover time of the smallest 
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turbulence eddies, as predicted by Kulsrud & Anderson (1992). During this kinematic 
growth period, the magnetic energy spectrum is peaked near k^p, and extends into the 
inertial range with a profile of ~ k 3 ^ 2 , providing seed magnetic field in the inertial range. 
Such seed magnetic field does grow, and the growth rate in the kinematic growth period 
is oc l/r max . Let range V include all the scales between k Vt D and k^o (see Table 1 for 
definitions of other terms). As the magnetic field within range V continues to grow, it starts 
to exert strong back reaction on the velocity field near scale k v> r> and suppresses it. This 
starts the dynamic growth period of the magnetic field. The eddy turnover time r max of the 
gradually suppressed velocity field increases, which reduces the growth rate, f3v ~ 1/T max , 
of the magnetic field within range V. Because (3 V is reduced, the back reaction of the 
magnetic field near the scales ~ k V:D on the velocity field is reduced, too. For the forced 
turbulence, line stretching in range X (the inertial range, see Table 1) is still in effect as 
the velocity field in this range is not completely suppressed by the growing magnetic field. 
Therefore, the magnetic field within the inertial range continues to grow at a rate smaller 
than l/r max but not smaller than l/r e ddy The growth rate of the total magnetic energy 
during the dynamic growth stage is a combination of the growth rates of E% and Ed, as 
shown as the e^ 1 stage of Figure 9. In Figure 11, we compare the growth history of E% and 
E v . Let £ = E x + E v , Ri = Ej/S and Rd = E v /£. Figure 11 shows that E x starts from 0, 
and continues to grow until it dominates Ex> after t ~ 40. In previous numerical simulations 
(Kinney et al. 2000), E% was not considered at all, and the growth of S is attributed only 
to E V - Figure 11 shows that E x is as important as E v in the amplification process of the 
magnetic field in turbulence, hence it must be included in numerical studies. 

We also studied a relevant process in Run D. The initial seed magnetic field is composed 
of two components: a large-scale magnetic field, B, and a small-scale magnetic field B fc=2 o 
that is initially concentrated at k = 20. The growth of the magnetic energy spectrum given 
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such combined initial conditions is shown in Figure 12. For this simulation run, we set the 
initial magnetic energy density of B equal to that of Bfe =2 o, so that we can compare the 
contributions of these two initial seed fields to the growth of the magnetic energy at each k. 
Figure 12 shows that both B and B fe=2 o contribute to the amplification of magnetic field at 
each scale. For a certain k within the inertial range, these two contributions to the magnetic 
energy density at this scale race against each other: the contribution from B initially grows 
as ~ t 2 , followed by an exponential growth ~ e */ Tedd f ; the contribution from B fc=2 o grows 
initially as ~ e *A"ma^ followed by another near-exponential growth with a smaller growth 
rate. The growth rate of magnetic energy density within the inertial range is greater than 
both of the growth rates in Run A and Run B, and this is because in Run D, both the 
large-scale and small-scale magnetic fields provide seed field within the inertial range. We 
believe this model can be applied to many real astrophysical systems, for example, in 
the regions where supernova remnants mix with ambient interstellar medium. Supernova 
remnants usually carry small-scale magnetic field, while the ambient interstellar medium 
can be threaded by large-scale magnetic field. Both the initial large-scale magnetic field 
in the interstellar medium and the small-scale seed magnetic field in supernova remnants 
will contribute to the growth of magnetic field at different scales. It is also possible that in 
the early evolution stages of galaxies, the seed magnetic field may have components of both 
large and small scales. As we have shown in above discussions, because the contributions of 
these two components to the amplification of magnetic field at different scales are different, 
we have to treat both of them in discussions of dynamo action in galaxies. 

Scale separation, which is discussed extensively in Brandenburg (2000), is not 
considered in this work. What we mean by "large scale" in this work is the scale k — > 0, 
while the "large scale" in Brandenburg (2000) is k ~ 1. That is, there are no scales between 
the large scale S and the forcing scale k ~ 1 in this work. We force the turbulence near 



-22 - 



k — 1, which is near the size of the simulation box, while the forcing scale of Brandenburg 
(2000) is k — 5. However, we did find the scale-by- scale growth of magnetic energy in our 
simulation. In the first two cases of our simulation (Runs A, B and C), in the kinematic 
development stage, magnetic energy near the kinetic dissipation scale is always stronger 
than the magnetic energy at different scales of the inertial range. When the magnetic field 
is amplified, the effect of back reaction first comes in near the kinetic dissipation scale and 
the growth of magnetic energy near that scale first slows down. After the slowing down 
of magnetic energy growth first appears near k v> D, the energy growth at the scale that is 
slightly larger than the dissipation scale starts to slow down, followed by the energy growth 
slowing down at even larger scales (or smaller fc's). Such slowing down, which starts from 
k Vi o, continues from the dissipation scale and moves from right to left in k— space, until 
the magnetic energy near the forcing scale stops growing. Then, the development enters 
fully saturated stage, with the magnetic energy at different scales fluctuating around their 
saturated values. Such successive slowing down of magnetic energy growth from the kinetic 
dissipation scale to the forcing scale is due to the effect of back reaction of the growing 
magnetic field on the velocity field, and such effect of back reaction comes in scale- by-scale, 
starting from k v> o- Note that this is different from the energy growth near scales k < 5 
in the work by Brandenburg (2000, Fig. 18) in the following sense: because the author 
is studying the dynamo a— and (3— effects, scale separation is necessary as the author 
compares the numerical results with the prediction of theoretical two-scale approach to the 
problem. In that work, the growth of magnetic field at k < 5 is mainly due to an a— effect, 
through which the helical velocity field aligns the magnetic field within the scales from 
k ~ 5 (i.e., the forcing scale) down to k V; n (i.e., the kinetic dissipation scale) in such way 
that the magnetic field at scale k ~ 1 is formed. In this work, on the other hand, the growth 
and the saturation of magnetic energy density at different scales is due to the interplay 
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between the magnetic field line stretching by the turbulent velocity field and the magnetic 
back reaction on the velocity field. 

6. Conclusions and future work 

We have carried out numerical simulations of 3D incompressible MHD with a magnetic 
Prandtl Pr ps 3. Both the initial large-scale seed magnetic field and the initial small-scale 
seed magnetic field can be amplified. For a large-scale seed magnetic field, the magnetic 
energy density grows as ~ t 2 for the first few turbulence eddy turnover times, followed 
by an exponential growth, of which the growth rate is oc l/r e dd y - For a seed magnetic 
field at an initial input small scale, during the kinematic development stage, magnetic 
energy can be transported to all scales larger and smaller than the initial input scale: 
near the magnetic dissipation scales magnetic energy is removed by magnetic resistivity, 
while from the outer scale of the turbulence to the initial input scale, the magnetic energy 
spectrum follows a profile of ~ k 3 ^ 2 . The measurement of the growth rate during this 
kinematic process confirms the theoretical prediction by Kulsrud & Anderson (1992) that 
the kinematic growth rate is oc l/r max where r max is the eddy turnover time of the smallest 
eddies of the turbulence. Entering the dynamic growth stage, the growth of magnetic 
field at small scales exerts a strong back reaction on the velocity field near the velocity 
dissipation scale. The suppression of the velocity field slows down the growth of magnetic 
field between the velocity dissipation scale and the magnetic dissipation scale. However, 
the magnetic field within the inertial range of the forced turbulence continues to grow. The 
magnetic field within the inertial range grows to a steady state that has a profile ~ k~ l 
between the forcing scale and the magnetic dissipation scale. The contribution to the total 
magnetic energy from the magnetic field within the inertial range dominates that from the 
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magnetic field between the velocity dissipation scale and the magnetic dissipation scales. 
For real astrophysical systems, the initial seed magnetic field may have both a large-scale 
component and a small-scale component, and they would both contribute to the growth of 
magnetic field at all turbulence scales. 

Our simulations may suffer from the relatively low resolution (64 3 ); therefore, the 
results may not completely applicable to very large magnetic Prandtl cases. The low 
resolution forces us to adopt the hyper- viscosity vj and the hyper-resistivity A7 in a few runs 
of our simulation. Although the introduction of hyper-viscosity and hyper-resistivity does 
not change the fundamental physics of our discussion given above, the simulation results, 
such as our estimates of the growth rates, can be slightly different from the results obtained 
from simulation runs with normal viscosity and normal resistivity, as we have discussed 
in Section 4. Nevertheless, our simulations show at least that the inertial range of the 
turbulence must be included in the discussion of how the growing small-scale magnetic field 
affects the growth of magnetic field at all scales of the turbulence. This is of paramount 
importance to the correct understanding of astrophysical dynamo processes. In future 
work, we expect to study the growth of magnetic field at all turbulence scales with larger 
numerical resolutions. 
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the organizers and all participants of the Astrophysical Turbulence Conference held at ITP 
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Table 1. Terms Used in This Paper and Their Physical Meanings 



What we call 


What we mean 


large scale 


scale ~ S, which — ► oo (or ks — » 0) in this work 


small scale 


a scale that is between k Vt o and kb t D 


outer scale of the turbulence 


scale ~ L 


forcing scale 


scale ~ L 


inertial range or range X 


scales between L and 1 Vi d 


range V 


scales between 1 Vi d and I^d {k v ,D < k < k^o) 


Ex 


magnetic energy density of range I 


Ex> 


magnetic energy density of range T> 


E v (b){t) 


history of kinetic (magnetic) energy density evolution 


E v{b) (k) 


kinetic (magnetic) energy spectrum 


Eb 


magnetic energy density of B 


eo(fc) 


magnetic energy density that is concentrated at scale k 
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Table 2. Measurements of various physical quantities in different simulation runs 



Run 


Eb 




e (k) 


V 






A 






Pr 


k v ,D 


h,D 


A 


5 x 10" 


-4 




hyper, v-j = 


5 x 10" 


-16 


hyper, A 7 = 


3 x 10" 


-20 


3 


13 


28 


B 






1 x 10~ 3 


hyper, v 7 = 


5 x 10" 


-16 


hyper, A 7 = 


3 x 10" 


-20 


3 


13 


28 


C 






1 x 10~ 5 


hyper, v 7 = 


5 x 10" 


-16 


hyper, A 7 = 


3 x 10" 


-20 


3 


13 


28 


D 


1 x 10" 


-3 


1 x 10^3 


hyper, v 7 = 


5 x 10" 


-16 


hyper, A 7 = 


3 x 10" 


-20 


3 


13 


28 


E 


5 x 10" 


-2 




normal, v — 


0.01 




normal, A = 


0.01 




1 


16 


16 


F 


5 x 10" 


-5 




normal, v — 


0.013 




normal, A = 


0.013 




1 


13 


13 
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Fig. 1. — Kinetic energy spectrum at R\ = 43 (circles) and Ra = 70 (stars) from pure 
hydrodynamic turbulence simulation with a resolution of 64 3 . We used normal dissipation, 
iA7 2 U, for these two runs. Also shown in this figure is the kinetic energy spectrum with 
a hyper- viscosity u 7 = 5.0 x 10~ 16 (triangles), from a 64 3 spectral simulation of pure 
hydrodynamic turbulence. 
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Fig. 2. — Kinetic energy spectrum E v {k) and magnetic energy spectrum Eb(k) at t — 0.3 for 
Run A, in which the initial large-scale seed magnetic field is B = Be y with a B = 0.0316. 
Both axes are plotted in logarithmic scales. For a kinetic energy spectrum of the form 
k~ p , the magnetic energy spectrum, which is generated from the interaction between the 
turbulence and B, follows k~ p+2 ~ k ' 77 . 



-31 - 





Fig. 3. — Pit), as defined in equation (13) in main text, at various times are plotted in the 
left panel of this figure. The averaged value (over wavenumber k) of (3(t) as a function of time 
is plotted in the right panel. In the left panel, the horizontal axis is plotted in logarithmic 
scale, while the vertical axis is plotted in linear scale. In the right panel, both axes are 
plotted in linear scales. Data shown here are from Run A, in which the initial large-scale 
seed magnetic field is B = Be y with a B = 0.0316. 
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Fig. 4. — Kinetic energy density E v = J2k E v (k) and magnetic energy density Eb = J2k Eb{k) 
as functions of time. Same conditions as Figure 3. The growth of magnetic energy density 
follows four stages. Stage 1 starts from the beginning of the simulation to t — 3, during which 
the magnetic energy grows as t 2 . During stage 2 that lasts from t — 4 to t — 16, it grows 
exponentially with growth rate f3 = 0.1. After the exponential growth stage, magnetic energy 
density growth enters a near saturation stage from t — 16 to t — 30, where the growth is 
further slowed down. For the saturation stage, the magnetic energy density fluctuates around 
0.33, while the kinetic energy density fluctuates around 0.4. 
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E b (k) at t=0.31 
E b (k) at t=0.94 
E b (k) at t=1 .57 
E b (k) at t=2.20 
E b (k) at t=2.83 
E b (k) at t=3.46 
E b (k) at t=4.71 
E b (k) at t=6.28 
E b (k) at t=7.85 
E b (k) at t=9.42 
E b (k) at t=21 .99 
E b (k) att=31.42 
E b (k) at t=40.84 
E v (k) at t=0.31 
E V (k) at t=4.71 



_ x _ E v (k) at t=9.42 

E^(k) at t=40.84 
^ E b (B y )=0.0005 
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Fig. 5. — Kinetic energy spectrum E^k) and magnetic energy spectrum E v (k) at different 
times. Data shown here are from Run A, in which the initial large-scale seed magnetic field 
is B = Be y with a B = 0.0316. The thick horizontal line denotes the energy density due 
to B. Growth of the magnetic field at all scales can be seen. The growing magnetic field 
suppresses the velocity field in the first three stages defined in Figure 4. 
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Fig. 6. — Isosurface of magnetic field magnitude at different times. The data are from 
Run A, in which the initial large-scale seed magnetic field is B = Be y with a B = 0.03f6. 
The isosurfaces are determined at 2.5 x mean{|£?|}, where mean{|_B|} = 0.23 for (a) and 
mean{|_B|} = 0.70 for (b). The emergence of magnetic structures of the size of the simulation 
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Fig. 7. — Magnetic energy spectrum E b (k) and kinetic energy spectrum E v {k) at various 
times. Data shown here are from Run B, in which an initial seed magnetic energy is 
concentrated at k = 20 with an e = 1 x 1CT 3 . Dashed lines are Ck 3 ^ 2 (C varies), which is 
reminiscent of the prediction by Kulsrud & Anderson(1992) for the magnetic spectrum that 
grows from a magnetic impulse. 
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Fig. 8. — Plots of the convective effect kCb(k) and the linear stretching effect kLb(k) at 
different time points. The horizontal axis is plotted in logarithmic scale, and the vertical 
axis is plotted in linear scale. Data shown here are from Run B, in which an initial seed 
magnetic energy is concentrated at k = 20 with an e = 1 x fCT 3 . Note the different scales 
of y-axes of different panels. 
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Fig. 9. — Temporal evolution of kinetic energy density and magnetic energy density. Data 
shown here are from Run B, in which an initial seed magnetic energy is concentrated at 
k = 20 with an e = 1 x 1(T 3 . 
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Fig. 10. — Magnetic energy density E b = Y,k Eb(k) as a function of time during the 
exponential growth stage of Run C, in which an initial seed magnetic energy is concentrated 
at k = 20 with an e = 1 x 10~ 5 . 7 b = 0.36. 7„ = 0.69 is the inverse of turnover time 
of the smallest eddies. 7 e // = 0.67 n as predicted by Chandran(1997) and Schekochihin & 
Kulsrud(2000). Redely = 0.18 is the inverse of turnover time of the largest eddies. 



-39- 




10 20 30 40 50 

t 



Fig. 11. — Comparison of the magnetic energy within the inertial range and the magnetic 
energy within the dissipation range. Note the linear-linear plot scales. Data shown here are 
from Run B, in which an initial seed magnetic energy is concentrated at k = 20 with an 
e = 1 x 1(T 3 . 
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Fig. 12. — Magnetic energy spectrum Eb(k) and kinetic energy spectrum E v (k) at different 
time points. Data shown here are from Run D, in which an initial seed magnetic field is 

2 

composed of a large-scale seed field, B = Be y with a B /2 = 1 x 1CT 3 , and a small-scale 
seed field that is concentrated at k = 20 with an e = 1 x 10~ 3 . Straight dashed lines are 
Ck 3 / 2 for various values of C. 



